On the ultra- high energy 
cosmic ray horizon 

D. Harari, S. Mollerach and E. Roulet 
CONICET, Centro Atomico Bariloche, 
Av. Bustillo 9500, Bariloche, 8400, Argentina 

February 4, 2008 

Abstract 

We compute the ultra-high energy cosmic ray horizon, i.e. the distance 
up to which cosmic ray sources may significantly contribute to the fluxes 
above a certain threshold on the observed energies. We obtain results 
both for proton and heavy nuclei sources. 

1 Introduction 

Soon after the discovery of the cosmic microwave background (CMB) radiation 
it was realized by Greisen, Zatsepin and Kuzmin [1] that the fluxes of cosmic 
ray (CR) protons with energies of order 10 20 eV and above would be strongly 
attenuated over distances of few tens of Mpc. This is due to the energy losses 
caused by the photo-pion production processes in the interactions of the protons 
with the CMB photons. Similarly, if CR sources accelerate heavy nuclei, these 
can photodisintegrate into lighter ones as they interact with CMB and infrared 
(IR) photons on their journey to us. In this way the fragments may arrive to the 
Earth with significantly smaller energies than the parent nuclei produced at the 
sources. Moreover, both protons and heavy nuclei can further loose energy by 
pair production processes, although due to the small inelasticities involved the 
typical attenuation length associated to e + e~ production at ultra-high energies 
is large (~ 1 Gpc for protons). 

Many works have studied these processes (see e.g. refs. [2, 3, 4, 5, 6, 7, 8, 
9, 10, 11] for the case of protons, and [12, 13, 14, 15, 16, 17, 18] for nuclei), 
and our purpose here is to make a detailed analysis to present the results in 
a way which can be useful for the study of correlations between the arrival 
directions of the highest energy cosmic rays (we will focus here on energies 
above 50 EeV, where 1 EeV = 10 18 eV) and candidate astronomical sources. 
Assuming that the sources are distributed uniformly and have similar absolute 
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cosmic ray luminosities, wc will obtain the distances within which the major part 
of the observed events should be produced (i.e. the horizon for the potential 
sources) as a function of the threshold adopted for the energies of the events. 



2 The horizon for protons 

The attenuation length A for the propagation of protons is just 

^ 1 = ~£ "^T = \l + Ke + K S C 1 ) 

receiving contributions from photo-pion (A 77r ), pair creation (A ee ) and redshift 
(X z ) losses. 

The attenuation length from photo-pion production by the protons in inter- 
action with the CMB photons is given by (we use hereafter natural units) 

A 77r = ^-J o de e V (e)a(e) 

= ~2^f I de'e'7 7 (e')a(e')ln[l-ex P (-e72 7 fcT)] ) (2) 

where 7 = E/m p is the Lorentz factor and the CMB photon density in the 
observer's rest frame is 

^ = t f3l 

de ^ 2 (exp(e/fcT)-l)' y> 

being characterized by a temperature T = 2.73°K. The cross section for pion 
production in terms of the photon energy e' in the proton's rest frame is a(e') 
and the inelasticity factor rj(e') is the average fraction of the proton energy lost 
in the process. We use for the product ar] a fit to the average values given in [8] , 
which includes both single pion production near the A resonance and multi-pion 
production at higher energies. 

This attenuation length, as well as that due to pair production obtained 
following [19, 20], are depicted in fig. 1 as a function of the proton energy. The 
redshift losses, characterized by X z ~ c/Ho ~ 4 Gpc, are of no relevance for 
the energies considered in the present work, for which the effects due to the 
interactions with the photon background are the dominant ones. 

A direct way to compute the CR horizon is to follow the trajectories of 
many individual protons with initial energies distributed according to a given 
input spectrum, here taken as a power law one, dN/dE oc E~ a (we will show 
results for a = 2.2 and 2.7). We then obtain the attenuation factor A(E t h,x), 
which is defined as the fraction of the particles which originally had an energy 
above a threshold value E t h that still have an energy above that threshold after 
traversing a distance x. Assuming that the sources are uniformly distributed and 
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Figure 1: Proton attenuation length vs. energy. 
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have similar intrinsic CR luminosities and spectra, one finds that the fraction of 
the events observed above a given energy threshold which originated in sources 
farther than a distance D is just 



F(E th ,D) = 



Jj di A(E th ,x) 
f™dx A{E th ,x)' 



(4) 



Since we will be interested in threshold energies above 50 EeV, cosmological 
effects due to the Universe's expansion or to source evolution are negligible. 
This implies that the effects of the inverse square distance reduction of the 
fluxes from each source and the increase in the number of sources with distance 
compensate each other, leaving just the simple integrals in eq. (4). 

The fraction F is depicted in fig. 2. We see that the horizon for protons, 
which may be taken as the distance for which this fraction reaches e.g. 10%, 
is relatively close on cosmological grounds for all the energies considered. For 
instance, for E t h = 80 EeV one has that 90% of the events should have been 
produced at distances not farther than <~ 90 Mpc 1 . The sensitivity to the 
assumed source spectral index a is not large, although as expected the horizon 
increases for harder spectra since above a given threshold the fraction of higher 
energy events, which are more penetrating, becomes larger in this case. Let us 
also mention that the pair production losses have an impact on the results only 
for E t h < 80 EeV, and they become indeed quite important for E t h < 60 EeV. 



1 It is not straight-forward to guess the results in fig. 2 from those in fig. 1, mainly because 
fig. 2 is in terms of the threshold energy measured on Earth, so that the CRs involved have 
energies above the thresholds and moreover their energies were even larger at the sources. 
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Fraction of events from x > D dN s /dE = k E" 



V. 'v, 1 1 

\ \ \ \\ 


i 


V\ \\ 
\ \ \ : '» '. 

\ \\ \0.1 




'. \ \\ V\ 
\\ \0.3\, 


a=2.7 (thick) 


:\ ; v. V- 

•.'\ ! V. V\ 

\\ v. \. 


<x=2.2 (thin) 


:V V* v-i : 




\\ V. X -L'- 




\s! •••X s ->.-. 
\Q / ••>,^ . 




i i i 





200 



150 



100 D[Mpc] 



50 



60 



80 



100 
E th [EeV] 



120 



140 



Figure 2: Considering the events arriving with energies above Eth, the con- 
tours indicate the fraction originating from distances larger than D, assuming 
a uniform distribution of proton sources with similar intensity. 



Note that since for a power law spectrum the number of events above a 
given threshold is proportional to E^ h a+ , one has simply that A(Eth,x) — 
(E t h/E' th (x)) a ~ , where E' th {x) is the initial energy that a CR proton should 
have in order to arrive with an average energy E t h after traversing a distance 
x. The energies E' th can be easily computed by solving the evolution equation, 
dE/dx = —E/X, backwards in time. We checked that this simpler procedure 
leads indeed to the same results as the Monte Carlo computations discussed 
before, but it has the drawback that it cannot be generalized to the case of 
heavy nuclei, since for these last it is not possible to know in advance the final 
average mass of the fragment after the nucleus has propagated a given distance. 

Let us mention that we are using the so-called continuous energy loss approx- 
imation to evolve the CR energies. This does not account for the fluctuations 
originating from the stochastic nature of the pion production processes. Al- 
though these fluctuations could be relevant for the study of detailed features in 
the differential CR spectrum (such as the 'bump' preceding the GZK cutoff), 
they have only a minor impact on the quantities we compute, which are based 
on the integral spectra above specified thresholds. 
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3 The horizon for nuclei 




The rate at which a nucleus of mass number A photo-disintegrates with the 
emission of i nucleons is given by 

1 f°° de An f 2 ~ ,e 

where now 7 = E/Am p . The relevant photon densities in this case are the 
CMB and the IR backgrounds. For this last we use the estimates obtained in 
[21]. The cross sections ua,i for the different processes in which a nucleus of 
mass number A emits i nucleons were parametrised in [12], and we also use 
the updated energy threshold values presented in [14]. To describe the average 
energy loss of a given nucleus it is convenient to introduce the effective rate [12] 

RA,eff = y^JRA,i- (6) 

i 

In terms of this rate one has 



The attenuation length for photodisintegration is then 

T d Adx A 
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The quantities A 7 d, together with the attenuation lengths to pair creation, are 
depicted in fig. 3 for different nuclei: 56 Fe, 28 Si and 4 He. It is apparent that pair 
creation losses become subdominant for all energies in the case of nuclei lighter 
than Si, while for instance in the case of Fe they are dominant for E ~ 100 EeV. 
On the other hand, for A > 20 one has that to a very good approximation 

\ ld (A,E)~\ ld (56,56E/A). (9) 

For light nuclei the attenuation lengths are in general somewhat larger than what 
would be implied by the above relation, due in particular to the less important 
role played by multi-nucleon emission in these cases and the narrower giant 
dipole resonances of lighter nuclei 2 . Anyway the general scaling of the results 
when different nuclei are considered can be understood from eq. (9). 



Fraction of events from x > D , Fe source , o!N s /dE = k E' 
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Figure 4: Considering the events arriving with energies above Eth, the con- 
tours indicate the fraction originating from distances larger than D, assuming 
a uniform distribution of Fe sources with similar intensities. 

In fig. 4 we depict the horizons for sources producing fluxes of Fe nuclei, 
obtained by following a large number of CRs with energies distributed according 
to the assumed spectra and obtaining the attenuation factors and associated 
fractions, as done before for the proton sources (fig. 2). For the evolution of 

2 Other peculiarities affect some nuclei, such as is the case for Beryllium, which disintegrates 
via Be —* 2 He + p. 
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Fraction of events from x > D , Si source , dNg/dE = k E 
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Figure 5: Considering the events arriving with energies above Eth, the con- 
tours indicate the fraction originating from distances larger than D, assuming 
a uniform distribution of Si sources with similar intensities. 

the mass number we use eq. (7) 3 , while for the energies we account also for 
pair creation losses. In the Fe source case the horizons result smaller than for 
proton sources and they become even smaller for lighter nuclei, for which the 
attenuation lengths become shorter at the energies we are considering (fig. 3). 
This is illustrated in fig. 5, corresponding to the case of 28 Si sources. It is clear 
that in a realistic situation the actual fractions will depend on the details of the 
non-homogeneous local source distribution and on the source intensities, but 
the limiting distances from which CRs can arrive above a given threshold are 
quite insensitive to these details. 

To better understand these results we show in figs. 6 the energies as a func- 
tion of distance for a set of 50 CRs with initial energies sampling an E~ 2 - 7 
spectrum for the three cases discussed (p, Fe and Si). The results in figs. 2, 4 
and 5 can then be interpreted naturally. 

An important difference between the proton and nuclei results is related to 
the fact that when a nuclear species is above the threshold of the giant dipole 
resonance for a typical CMB photon, also the fragments from the photodisinte- 
gration will be above threshold (neglecting the subdominant energy losses due 
to pair creation) and hence the disintegration continues efficiently until the nu- 

3 For a given mass number wc take the charge as the smaller one corresponding to a stable 
isotope. 
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Figure 6: Energies as a function of distance to the source for a sample of initial 
energies. Top figure is for protons and bottom ones for Fe and Si sources (note 
the different horizontal ranges). 



cleus is completely disintegrated. On the contrary, nucleons which are originally 
above the threshold for photo-pion production and loose energy by pion emis- 
sion may soon end up below threshold and hence continue their travel without 
further losses due to pion production. This leads to the well known pile-up in 
the proton spectrum, which is well seen in the top panel of fig. (6), a feature 
which is much less pronounced in the case of nuclei. This effect is also respon- 
sible for the smaller sensitivity of the results in figs. 4 and 5 to the adopted 
spectral index. 

It is interesting to note that even if these scenarios adopt a single nuclear 
mass at the sources, the distribution of the masses of the CRs arriving to the 
Earth is quite wide, with contributions from many different nuclear species. This 
is exemplified in fig. 7, which displays the relative distribution of mass fragments 
arriving to the Earth above a threshold energy of 60 EeV for scenarios with Fe 
(dashed histogram) and Si (solid histogram) sources, with spectral index a = 2.7 
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(the results are not very sensitive to the value of a adopted). 

The flux of secondary nuclcons produced in the various photodisintcgration 
processes is generally quite suppressed, since these nucleons will be emitted 
with energies E/A if the CR started with energy E and mass A. Hence, the 
secondary nucleons can contribute non-negligibly to the CR fluxes only in the 
case of very hard spectra (a ~ 2) and high source cutoffs (E max ^> AE^). For 
simplicity, we have just assumed in our computation that the maximum energy 
attained in the sources is E max = AE th (which is bigger than 10 21 eV in the 
examples considered), and in this case the secondary nucleons are always below 
threshold. 
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Figure 7: Distribution of masses of the fragments arriving to the Earth with 
E > 60 EeV in a scenario with uniformly distributed Fe (dashed) and Si (solid) 
sources with spectral slope a = 2.7 (arbitrary units). 

It is also worth mentioning that, as is seen from fig. 3, heavy nuclei produced 
with energies E > A x 10 EeV have attenuation lengths of only a few Mpc. 
Hence, these nuclei will be essentially completely disintegrated after traveling 
~ 10 Mpc. Since each nucleus of initial mass A and energy E will lead to A 
nucleons with energies E/A, the number of nucleons above an energy threshold 
E will be in this case N(> E ) ~ AN A (> AE ), with N A (> E) cx E 1 - 01 being 
the original integral source spectra (assuming here that the power law flux has 
no upper cutoff). We hence learn that the sources of nuclei farther than 10 Mpc 
will produce a flux of secondary nucleons above 10 EeV similar to the source flux 
but suppressed by a factor A 2 ~ a , and extending up to energies E max /A. The 
attenuation of these secondary nucleons caused by photo-pion and pair creation 
losses may then be studied as done in section 2. 

Let us finally mention that both pair production losses and photo-disintegrations 
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with the IR background have a significant impact, especially for the lower thresh- 
olds considered, on the results for Fe nuclei (fig. 4), while in the case of Si sources 
(fig. 5) their effect is less important. 

4 Conclusions 

As a summary, we have shown that the sources of the vast majority of the 
CR events observed with energies above 50 EeV should be at relatively close 
distances (under the plausible assumption that CRs consist of ordinary nucleons 
or nuclei), as is illustrated in a quantitative way in figs. 2, 4 and 5 4 . These 
results should be relevant in order to restrict the distances to the candidate 
astrophysical objects that could act as sources for the observed ultra-high energy 
cosmic rays. 
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